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The surprising recent discoveries of quasicrystals and their approximants in soft matter systems 
poses the intriguing possibility that these structures can be realized in a broad range of nano- 
and micro-scale assembhes. It has been theorized that soft matter quasicrystals and approximants 
are largely entropically stabilized, but the thermodynamic mechanism underlying their formation 
remains elusive. Here, we use computer simulation and free energy calculations to demonstrate 
a simple design heuristic for assembhng quasicrystals and approximants in soft matter systems. 
Our study builds on previous simulation studies of the self-assembly of dodecagonal quasicrystals 
and approximants in minimal systems of spherical particles with complex, highly-specific interac- 
tion potentials. We demonstrate an alternative entropy-based approach for assembling dodecagonal 
quasicrystals and approximants based solely on particle functionalization and shape, thereby re- 
casting the interaction-potential-based assembly strategy in terms of simpler-to- achieve bonded and 
excluded-volume interactions. Here, spherical building blocks are functionalized with mobile surface 
entities to encourage the formation of structures with low surface contact area, including non-close- 
packed and polytetrahedral structures. The building blocks also possess shape polydispersity, where 
a subset of the building blocks deviate from the ideal spherical shape, discouraging the formation 
of close-packed crystals. We show that three different model systems with both of these features 
- mobile surface entities and shape polydispersity - consistently assemble quasicrystals and/or ap- 
proximants. We argue that this design strategy can be widely exploited to assemble quasicrystals 
and approximants on the nano- and micro- scales. In addition, our results further elucidate the 
formation of soft matter quasicrystals in experiment. 



Until fairly recently, quasicrystals and their approx- 
imants have been observed only in atomistic systems. 
Over the past decade, there have been sporadic reports 
of quasicrystals and approximants in nanometer and 
micron-scale systems. Examples include holographically- 
trapped T and laser-field-induced [2 [3] quasicrystalline 
materials made of micron-sized spheres, self- assembled 
quasicrystals and approximants formed by spherical 
dendrimer micelles [H |3], phase-separated star-triblock 
copolymers [5], binary nanoparticle superlattices [7], 
spherical micelles of phase-separated block copolymers [3 
[9], and simulations of hard tetrahedra [10]. These re- 
ports pose an intriguing possibility that these structures 
might be assembled in a broad range of systems. In one 
such system, spherical dendrimeric micelles functional- 
ized with alkyl tails form a dodecagonal (12-fold) qua- 
sicrystal (DQC), as well as other non-close packed struc- 
tures such as the body-centered cubic (bcc) and A15 
crystals In similar systems, various types of block 

copolymer micelles arrange into quasicrystals with 12- 
fold, and possibly 18-fold, symmetry [S":, as well as vari- 
ous periodic approximants 8 . 

The dendrimer and block copolymer micelle systems 
in particular all share an important common feature: 
their constituent micelles exhibit a soft "squishy corona" 
in which terminal groups avoid each other to mini- 
mize steric interactions. It has been predicted that this 
mechanism causes the system to adopt structures that 



minimize surface contact area between neighboring mi- 
celles [12l [13] . The structure that minimizes surface con- 
tact area, known as the Weaire-Phelan or A15 struc- 
ture [13], is structurally similar to a DQC, but, since 
DQCs do not minimize surface contact area, other fac- 
tors must contribute to their stability. It has been sug- 
gested that entanglement of terminal groups may give 
rise to three-body entropic effects that favor quasicrys- 
tals in systems of monodisperse micelles [THl [E] . In all 
these micellar systems, entropic effects appear to play a 
predominant role in stabilizing the quasicrystals and ap- 
proximants, potentially distinguishing them from many 
of their atomistic counterparts in which strong attractive 
interactions are present. 

Computer simulation studies of self-assembly have 
demonstrated that quasicrystals can be assembled by an 
inverse-design mechanism. In particular, pair potentials 
can be designed to make close-packing unfavorable, caus- 
ing such systems to instead form quasicrystals and ap- 
proximants |17H19j . These complex interaction poten- 
tials have yet to be realized in experimental systems on 
the micro- or nano-scale, but we propose that a similar 
effect can be achieved via shape polydispersity, where 
a subset of the micelles deviate from the ideal spheri- 
cal shape. Shape polydispersity arises naturally in many 
micelle-forming systems, and, in general, particle shape 
is a tunable parameter in many micro- and nano-scale 
systems [20) . 



In this article, we introduce a design strategy based on 
the ideas described above to direct the self-assembly of 
three-dimensional DQCs and/or their periodic approxi- 
mants in systems of (approximately) spherical micelles 
or similarly shaped particles. We study different types of 
nano/microscale building blocks with features that pro- 
mote structures with low surface contact area and sup- 
press close-packing. Structures with low surface contact 
area are promoted by functionalizing spherical building 
blocks with mobile entities connected to their surface, 
similar to functionalized spherical dendrimers [5]. Close- 
packing is suppressed by incorporating shape polydisper- 
sity into the system in the form of particle asphericity. 
Both features are relatively common aspects of soft mat- 
ter and related systems and should be achievable ex- 
perimentally; a schematic of our strategy is shown in 
Fig. la. Applying this strategy in computer simula- 
tions, we show three key results. (1) We verify the the- 
oretical predictions that interactions between terminal 
coatings can drive the system to form surface-area min- 
imizing structures [T^l US]- (2) We demonstrate that 
shape polydispersity can be used to suppress the forma- 
tion of close-packed structures. (3) We show that three 
different simulated micellar systems that possess both 
of these characteristics reproducibly form DQCs and/or 
approximants. These models - a simplified model of a 
spherical micelle and two micelle-forming systems com- 
posed of tethered nanosphere building blocks piH^ - 
represent the only simulated micellar systems currently 
known to form 3d quasicrystals or approximants through 
self-assembly. Because the models are closely related to 
experimental systems known to form DQCs and/or ap- 
proximants [H IHl [H] , our results may provide pertinent 
insight regarding their formation. In the future, the as- 
sembly strategy that we employ may serve as a heuristic 
for expanding the range of systems that assemble DQCs 
and approximants. 



DODECAGONAL QUASICRYSTALS AND 
APPROXIMANTS 



We first introduce definitions and terminology that will 
facilitate our discussions in subsequent sections. A crys- 
tal is defined as a structure with long-range positional or- 
der, as identified, for example, by the presence of Bragg 
peaks in the diffraction pattern [26]. A quasicrystal is a 
quasi-periodic crystal; that is, a crystal that lacks period- 
icity [57], but still exhibits diffraction peaks. Quasicrys- 
tals sometimes (but need not) exhibit rotational sym- 
metries that are incompatible with periodicity. Several 
types of quasicrystals have been observed in experiment, 
but in this article we focus on DQCs in particular be- 
cause those are to date the most commonly reported type 
of quasicrystal in soft matter systems. DQCs are char- 
acterized by their long-range dodecagonal (i.e, 12-fold) 
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FIG. 1: Assembly strategy and structure of the dodecago- 
nal quasicrystals and approximants. (a) Schematic of the 
proposed two-part strategy that uses functionalization and 
shape to form DQCs. Particle functionalization (left) pro- 
motes the formation of structures with low surface contact 
area and asphericity (right) inhibits the formation of close- 
packed structures. Particles colored red in the asphericity 
schematic (right) are meant to highlight where the crystal is 
disrupted by the presence of aspherical particles (blue), (b) 
Valid tiles for the DQC. The DQC and approximants can be 
described as a periodic stacking of plane-filling arrangements 
of tiles in the z direction (out of the page). The gray particles 
at the nodes of the tiles form layers at z=l/4 and z=3/4 and 
sit at the centers of 12-member rings. The yellow particles and 
red particles form layers at z=0 and z=l/2 respectively. In the 
DQC, the gray particles form a dodecagonal layer with 12-fold 
symmetry, and the yellow and red particles form hexagonal 
layers rotated by 30 degrees to obtain 12-fold symmetry, (c) 
Three common DQC approximants. (d) A higher-order ap- 
proximant generated through the inflation method (see text) . 
(e) A representative DQC random tiling of squares, triangles, 
rhombs and shields, adapted from Ref. [25] 



rotational symmetry. 

DQCs are polytetrahedral structures [55] of the Frank- 
Kasper (FK) type [53]. For the class of FK structures 
considered here, ordered structures are distinguished by 
their "tiling" pattern, constructed by connecting the cen- 
ters of neighboring 12-member rings of particles (see pan- 
els b-e of Fig. 1). The structures are layered and, whether 
periodic or aperiodic in the plane, they repeat periodi- 
cally in the direction orthogonal to the plane (into to the 
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page in Fig. 1). There are five valid tiles that can be 
arranged to form structures with complete 12-member 
rings without disorder. These tiles take the shape of a 
square, a triangle, a rhomb, a shield, and an asymmetric 
hexagon jEj, and are illustrated in Fig. lb. Periodic 
arrangements of these tiles result in periodic crystals, 
sometimes known as "approximants," that are indistin- 
guishable from DQCs locally [30]. Three common ap- 
proximants, known as the A15, Z, and sigma structures, 
are shown in Fig. Ic. Increasingly complex approxi- 
mants, such as the structure depicted in Fig. Id, can be 
constructed by inflation, whereby tiles are sequentially 
replaced with smaller sub-tiles (STJ [32] . 

In addition to periodic arrangements, non-periodic ar- 
rangements of tiles that fill the plane can also be con- 
structed, resulting in quasicrystals. Various methods can 
be used to construct the tilings; methods such as in- 
fiation [32], projection ^T\, or matching rules [521 131] 
produce deterministic quasicrystals, whereas random 
tilings [3S] give rise to a range of similar quasicrystals 
with the tiles reshuffled locally, characterized by local 
phason fluctuations. Imperfect quasicrystals of either 
type may also exhibit global phason strain whereby par- 
ticular tiles or orientations of tiles occur more or less 
frequently that in the ideal case, giving rise to shifts 
and broadening of the diffraction peaks [36] . Determin- 
istic quasicrystals are thought to be energetically stabi- 
lized, whereas random tiling quasicrystals are thought 
to be entropically stabilized [3S]- Fig- le shows a typi- 
cal random-tiling DQC [55] that we envision might form 
in soft-matter systems, which are often stabilized by en- 
tropy. The structure is composed mostly of squares and 
triangles, and is locally similar to the sigma approximant. 
The sigma approximant is the thermodynamically sta- 
ble state for many systems that form DQCs, and the 
two structures often arise in nearby regions of parame- 
ter space [11 [5l [7] . The experimental protocol may dic- 
tate whether a metastable DQC or a stable sigma ap- 
proximant is obtained. In the case of the simulations 
we perform on model micelles, we are limited to rela- 
tively small, finite size simulations, as discussed subse- 
quently. As such our systems are typically too small to 
unambiguously distinguish between quasicrystals and ap- 
proximants, or identify phason strain. With this caveat 
in mind, we refer to our assembled structures as qua- 
sicrystals if they are composed of valid tiles for the DQC, 
exhibit strong peaks in the diffraction pattern, and are 
not periodic (aside from the trivial periodicity imposed 
by the periodic boundary conditions on the scale of the 
sample) . 

SIMULATION RESULTS 

We begin by performing molecular dynamics simula- 
tions [3S] of a simplified model of a spherical micelle 



(MSM) that considers only excluded volume interactions 
between terminal groups on the micelle surface (Fig. 2a) . 
Unlike the truly minimal "fuzzy sphere" micelle model 
of Ref. [T2] that treats inter-micelle interactions with an 
effective pair potential, our model treats these excluded 
volume interactions explicitly through mobile spheres at- 
tached to the micelle surface. This allows us to (1) study 
the self-assembly of the micelles and (2) directly mea- 
sure the relative effect of entropy and energy in driving 
the stabilization of assembled phases. The MSM consists 
of a non-interacting rigid scaffolding with 42 points on 
the surface of a sphere, given by the vertex points of a 
2-frequency icosahedral geodesic with diameter — 5.27a. 
With this diameter, the average spacing between surface 
points is 1.5(7. Each surface point anchors a small spher- 
ical particle with diameter a. The particles and surface 
points are attached by harmonic springs of stiffness k 
that control the degree of surface particle mobility. Sur- 
face particle mobility increases as k decreases, creating a 
larger, "squishier" outer corona. Decreasing k can also 
be interpreted as increasing the radius of gyration of the 
surface coating, if we consider the spheres to be dumb- 
bell polymers anchored to the surface [35] . Excluded vol- 
ume interactions between the surface spheres are mod- 
eled by the purely-repulsive Weeks-Chandler-Andersen 
(WCA) potential [ID] (see Materials). Roughly speaking, 
the MSM can represent many different nanoscopic ob- 
jects, including core-satellite nanoparticlcs [¥TH15] . where 
nanospheres are functionalized with an outer coating of 
smaller nanospheres; spherical micelles composed of den- 
drimers [5] [T^l 113] where the outermost layer of the den- 
drimer "tree" is functionalized with oligomers or poly- 
mers; spherical block copolymer micelles [5] [3] that pos- 
sess an outer corona of polymers; or spherical micelles 
created from amphiphilic tethered nanoparticlcs [211123] . 
as we discuss later. 

In the absence of shape polydispersity, the MSMs tend 
to form close-packed (face centered cubic, fee, or hexag- 
onally close packed, hep) arrangements for k > 5 (lower 
surface particle mobility) and bcc structures for fc < 5 
(higher surface particle mobility); structures are identi- 
fied using the algorithms described in Ref. [H] . These re- 
sults support the conjecture that increasing surface par- 
ticle mobility drives the system towards structures with 
lower surface contact area (such as bcc), as we discuss 
in detail in the following section. At these statepoints, 
sphere packing constraints favor the bcc structure over 
the surface-contact-area-minimizing A15 structure. A 
bcc-ordered structure of 60 MSMs is shown in Fig. 2b 
for k=5. 

We find a more dramatic change in the structural ar- 
rangement of the MSMs when shape polydispersity is 
incorporated into the system in the form of aspherical 
"dimer" micelles (see Fig. 2c). We allow dimers to form 
in an unbiased manner by exploiting the fact that at low 
fc, surface particles are only loosely bound to the sur- 
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FIG. 2: MSM systems, (a) MSM monomer (white) extracted 
from a simulation, (b) 60 MSMs at fc = 5 with /dimcr = 
(no dimers) in a bcc configuration, (c) MSM dimer (cyan) 
extracted from a simulation, (d) 60 MSMs with fc = 5 and 
/dimer ~ 0.24 lu a slgma structure, (e-g) Systems with 360 
MSMs and (e) k =4 and /dimcr ~ 0.39 , (f) k =4.75 and 
/dimer ~ 0.37 , aud (g) fc =5 and /dimcr ~ 0.36. In all cases, 
we plot time-averaged density isosurfaces of the centers of 
mass of the micelles/dimers, rather than the micelles them- 
selves, to remove thermal noise and produce a clearer image 
(see Materials). Systems are viewed along the pseudo-12-fold- 
symmetry axis, as calculated using the diffraction pattern, 
shown to the right of each panel. In all cases, systems are 
colored-coded with monomer centroids shown in white/grey 
and dimer centroids shown in cyan. Note, (e) appears as a 
parallelogram due to the projection; all simulation boxes are 
square or rectangular cuboids. 



face sites on the scaffold, allowing the MSMs to overlap; 
some of the MSMs become locked together into dimers 
when k is increased. By slowly increasing from a highly- 
disordered state at fc = 2, we create systems with dimer 
fraction in the range 0.20 < /dimer < 0.40, consisting 
of dimers with an average aspect ratio of ~1.45:1. This 
procedure roughly mimics the process by which micelles 
are formed in amphiphilic soft matter systems, such as 
the tethered nanoparticle models that we discuss later. 
In such systems, spherical micelles assemble from a disor- 
dered mixture of individual building blocks as the system 
temperature is reduced [32 03] (see Materials) . In the 
MSM system, increasing fc has a similar effect to decreas- 
ing the temperature. 

We find that systems with a mixture of spherical and 
dimer MSMs consistently form FK structures Fig. 
2d shows a typical sigma approximant formed by 60 
MSMs at fc = 5 with /dimcr = 0.24; sigma structures were 



reproducibly observed in over 25 independent simulations 
where fc was slowly increased from 2 to 5. This approxi- 
mant closely matches the expected result for 60 particles 
interacting via the Dzugutov or Lennard- Jones-Gauss 
pair potentials at densities that yield DQCs for larger 
systems. The formation of the sigma structure is also 
consistent with the observed experimental behavior of 
spherical dendrimer [4] and block copolymer micelles [8] . 
Three representative independent simulations, each com- 
posed of 360 MSMs in rectangular boxes with aspect ratio 
1.28:1.28:1.00, are presented in Figs. 2e-g. Figs. 2e,f, and 
g show systems at fc =4, 4.75, and 5, with /dimcr =0.39, 
0.37, and 0.36, respectively. In all cases, we observe 
finite-size DQCs that exhibit long-range rotational or- 
der of the MSM center-of-mass but no periodicity aside 
from the trivial periodicity imposed by the boundary con- 
ditions. Our simulations are limited to smaller system 
sizes than typical point-particle models [ITl [19] because 
we must resolve timescales corresponding to the micro- 
scopic motions of the surface particles that comprise the 
MSM, rather than the MSM centroid. Nevertheless, the 
finite structures depicted in Fig. 2e-g exhibit local in- 
dicators of DQC ordering. The systems form unique 
tilings with different configurations rather than any par- 
ticular approximant. The systems also contain the en- 
tire range of valid tiles, rather than containing squares 
and triangles exclusively like the sigma phase, which of- 
ten competes with DQCs for stability. Since DQCs grow 
more easily than approximants [37], it is possible that 
the DQC-like tilings are thermodynamically metastable 
relative to a stable approximant. The structures do not 
rearrange or undergo phason flips after solidification dur- 
ing the timescale of our simulations. 

We can further test our proposed strategy in systems 
where we do not have explicit control over surface par- 
ticle mobility or shape polydispersity, but where these 
two key features instead emerge naturally as a result 
of phase separation. We consider two model tethered 
nanosphere (TNS) systems, mono-TNS [22 [12] and di- 
TNS [13 [H [M], both of which form roughly spher- 
ical micelles with mobile surface entities. Schematics 
of the building blocks are shown in Figs. 3a, e respec- 
tively, and the micelles they form are shown in Figs. 3d,g 
respectively. The mono-TNS micelles have an outer 
shell of mobile nanospheres that closely match the MSM 
model, while the di-TNS micelles have a shell of short 
polymers, more closely resembling the spherical micelles 
formed by block copolymers [8j |46l [47] and functional- 
ized dendrimers[51 [T^ [T5]. These models are computa- 
tionally expensive, and thus only relatively small systems 
in terms of the number of micelles are explored. Fig. 3b 
depicts density isosurfaces [48] of the aggregating teth- 
ers for a system of 2500 mono-TNS building blocks that 
assemble into ~60 spherical micelles arranged in a sigma 
approximant. Fig. 3c depicts isosurfaces for a system of 
5000 mono-TNS that self-assemble into ^-^120 spherical 
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FIG. 3: TNS systems, (a) Schematic of a mono-TNS building 
block, where the 8 tether beads (blue) of size a aggregate and 
self-assemble spherical micelles with a soft core surrounded by 
relatively hard "satellite" nanoparticles (white) of size 2.5ct 
that act as mobile surface entities, (b) A simulation snapshot 
of ~60 micelles formed by mono-TNS that arrange into a 
sigma approximant, and (c) ~120 mono-TNS micelles that 
form a DQC-like structure; for both systems (j) ~ 0.275 and 
T=l.l. (d) Histogram of asphericity, as, of the mono-TNS 
micelles in the sigma phase, (e) Schematic of the di-TNS 
building block, where the 4 beads in the tether (green) each 
of size a aggregate and nanoparticles (white) of size 2a are 
also attractive; 4 bead tethers (purple) of bead size a that do 
not aggregate coat the outside of the micelle, (f) ~60 di-TNS 
micelles arranged in a sigma approximant at — 0.2 and T 
— 1.2. (g) Representative di-TNS micelles with different as- 
In all cases, for clarity, we show density isosurfaces of the 
aggregating polymer tethers (i.e., the micelle core). 



micelles arranged in a FK structure containing squares, 
triangles, shields, and rhombs. The increasing complex- 
ity of the tiling arrangement with system size indicates 
that the TNS system may form a higher-order approx- 
imant or a DQC in the infinite limit. The mono-TNS 
micelles naturally exhibit shape polydispersity. Fig. 3d 
shows a histogram of the asphericity, , computed from 
the principle radii of gyration |45j of the micelles, with 
representative micelles at various values of Os depicted in 
the inset. For reference, as — corresponds to an ideal 
sphere and Ug = 0.02 corresponds to the MSM dimer 
with aspect ratio 1.45:1 shown in Fig. 2c. Fig. 3f shows a 
sigma structure formed from 2000 di-TNS building blocks 
that self-assemble into ^60 micelles. The distribution of 
Us for the di-TNS (plotted in the Fig. 6) is similar to 
that for the mono-TNS system. Two representative di- 
TNS micelles at low and high are depicted in Fig. 3g. 
Overall, FK structures self-assembled from TNS build- 



ing blocks were reproducibly observed in 20 independent 
simulations. Whether these systems form DQCs in the 
infinite limit remains an open question that should be 
explored in the future as computational power increases. 

FREE-ENERGY CALCULATIONS 

Having explored the self-assembly of the three micelle 
models, we now perform free-energy calculations to inves- 
tigate the thermodynamic basis underlying both aspects 
of our strategy for DQC-like structure stabilization. The 
first aspect, the functionalization of particles with mo- 
bile surface entities, is inspired by the observation that 
soft-matter systems with relatively soft inter-micelle in- 
teractions often form non-close packed structures, as de- 
scribed in the Introduction. For example, spherical den- 
drimeric micelles functionalized with alkyl tails to create 
a "squishy corona" are known to form non-close packed 
structures such as the bcc and A15 crystals [TT]. Ziherl 
and Kamien proposed that the formation of the bcc and 
A15 structures is related to the Kelvin problem, which 
involves finding the space- filling arrangement of cells that 
minimizes surface contact area |12| I13j. In this picture, 
the dendrimeric micelles adopt structures with low sur- 
face contact area in order to reduce steric interactions 
between terminal polymer groups. The bcc and A15 crys- 
tals both exhibit low surface contact area, with A15 rep- 
resenting the current best-known solution to the Kelvin 
problem [T3]. It has been suggested [S] that this mech- 
anism may also stabilize the dendrimer DQC observed 
in experiments [S]. However, since minimizing surface 
area alone favors the A15 structure rather than the DQC, 
other factors must be important as well. 

We calculate the Helmholtz free energy, F |SD] , as a 
function of the surface particle mobility k for a system of 
monodisperse MSMs (i.e., without dimers); see the Ma- 
terials for more information. The value of F in Fig. 4a is 
shown relative to the value for the hep crystal, taken as 
a convenient reference state. Fig. 4a illustrates that as 
k decreases (i.e., surface particle mobility increases), F 
decreases more rapidly for the A15, dod, and bcc struc- 
tures than for the fee and hep structures. Here, the value 
for the "dod" curve is the average of the sigma phase and 
several higher-order square-triangle approximants to the 
DQC all of which have nearly identical free ener- 
gies. For low fc, bcc appears to be the stable state, con- 
sistent with our MSM simulation results. For very low k 
[k < 3) the system becomes disordered. The change in 
F as a function of k is the strongest for the A15 struc- 
ture, which minimizes surface contact area, followed by 
the dod and bcc structures, respectively. We note that 
the dod structure has a lower free energy than the A15 
structure over the entire range; however, at sufficiently 
low fc, the difference in free energy between bcc, A15, 
and dod is indistinguishable. The change in F with k 
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FIG. 4: (a) Helmholtz free energy per micelle, F, as a function 
of surface particle mobility (controlled by the spring stiffness 
k) for monodisperse MSMs. The inset shows the potential 
energy per micelle U. (b) _F as a function of /dimor for the 
Weeks-Chandler- Andersen sphere/dimer mixture. For parts 
(a) and (b), the energies and free energies are reported with 
respect to the hep crystal for convenience. For all datapoints, 
error bars are smaller than the data labels, (c) Depiction 
of three different FK polyhedra. (d) Probability of observing 
dimers at the center of Z12, Z14, and Z15 structures in the dod 
phase as a function of number density, p. (e) Fraction of Z12, 
Z14, and Z15 local structures in A15, sigma and Z structures. 
Note, (d-e) are color coded following the convention in (c). 



is entropically driven; the difference in average potential 
energy (U) changes little, and does not decrease with F 
(Fig. 4a, inset). This serves as a direct verification of 
the predictions of Ziherl and Kamien [12, J^, . Note that 
the Z structure (Fig. Ic) is omitted as it is not stable 
in the parameter range under consideration. While the 
trends in entropy are as we expect, we find surface par- 
ticle mobility alone is not sufficient to stabilize DQCs or 
approximants for the statepoints and model under con- 
sideration. Thus, as our self-assembly simulations pre- 
viously showed, a second mechanism is needed to form 
DQC structures for this model. 

The thermodynamic basis underlying the second as- 
pect of our strategy - shape polydispersity - can be 
understood in the context of previous studies of both 
quasicrystal formation and sphere packing. Systems of 



particles with short-ranged, spherically-symmetric inter- 
action potentials, such as hard spheres or particles with 
short-ranged van der Waals interactions, modeled by the 
Lennard- Jones (LJ) potential, tend to form close-packed 
crystals in the solid phase, e.g. f and/or hep. Al- 
though these systems tend to locally favor polytetrahe- 
dral structures |51) , close-packed structures maximize the 
overall packing density and hence maximize the entropy, 
and also often exhibit low potential energy. Special- 
ized interparticle potentials, such as the Dzugutov [17] 
and Lennard- Jones-Gauss [T^ potentials, have been con- 
trived with features that help drive systems away from 
close-packed structures. Like the standard LJ potential, 
the Dzugutov and Lennard- Jones-Gauss potentials have 
an attractive well that encourages local polytetrahedral 
ordering. However, these specialized potentials include 
an additional relative energy penalty for adopting the 
characteristic interatomic spacings associated with close 
packing, ultimately driving the system to form alterna- 
tive structures, such as bcc crystals, as well as DQCs and 
their approximants under certain conditions [THJ [T5] . We 
propose, as our previous MSM simulations show, that 
shape polydispersity can have a similar effect, driving the 
system away from close-packing. However, in contrast to 
the energetic repulsion of the Dzugutov potential, the 
destabilizing effect, in this case, is entropic. 

To explicitly quantify the effect of shape polydispersity, 
we perform free energy calculations |52H56j for binary 
mixtures of soft spheres and short, pill-shaped dimers, 
with particle interactions modeled by the WCA poten- 
tial (see Materials). The dimers are modeled by a rigid 
body of length 1.5a consisting of two overlapping soft 
spheres 0.5a apart (see Fig. 4b), resulting in an aspect 
ratio of 1.5:1, similar to the aspect ratio observed in the 
simulation of MSMs. Fig. 4b shows the Helmholtz free 
energy, F, as a function of the dimer fraction, /dimor, 
for several structures at a representative state point with 
number density p = 0.9 and T — 0.25. The free en- 
ergy is computed based on the standard Einstein crystal 
thermodynamic integration (TI) method [SH [5S], with 
an additional alchemica l[56j TI step to compute the free 
energy required to transform a given fraction of spheres 
into dimers (see Materials). As /dimcr increases, the A15 
and dod structures become increasingly stable relative 
to close packed crystals, and, to a lesser extent, the bcc 
crystal. We note that the dod phase has a lower value 
of F than the A15 structure for all statepoints, although 
the difference becomes minimal for high dimer fraction. 

This difference in stability between the FK phases (A15 
and dod) and standard crystals can be traced to the ten- 
dency for dimers to adopt larger, more aspherical neigh- 
bor shells, which are present in FK structures but not fee, 
hep or bcc crystals. The first neighbor shells of particles 
in FK structures form different types of polyhedra, which 
may be icosahedral (coordination number 12), or take on 
higher coordination numbers Z, such as Z13, Z14 or Z15 
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depicted in Fig. 4c. In Fig. 4d, we plot the probability of 
observing dimers in Z12, Z14, and Z15 configurations for 
the dod structure where we fix particle centroids but al- 
low dimers to rotate and swap positions with monomers. 
We observe that dimers strongly favor Z15 coordination 
shells as these are the largest and thus most accommo- 
dating. Dimers sit in ZI4 arrangements as a second best 
option and almost never occupy Z12 structures which are 
the smallest. We can gain additional insight by examin- 
ing the relative fraction of Z12, Z14, and Z15 within the 
three approximant structures, as shown in Fig. 4e. Al- 
though the free energy of the A15 and dod phases are 
similar, the A15 phase does not possess any Z15 arrange- 
ments, whereas the dod phase has an appreciable fraction 
(^0.13). This difference may account for the widespread 
formation of dod rather than the A15 structures in our 
three simulation models. We note that while the Z phase 
has the largest fraction of Z15 coordinations, it also pos- 
sesses the largest fraction of the less favorable Z12 co- 
ordinations, which may partially account for its relative 
instability for this density and dimer size. 

We observe that for /dimor > 0.4, A15 and dod are 
more stable than fee, hep and bcc crystals. This implies 
that mixtures of spherical and pill-shaped colloids might 
produce DQCs or approximants. However, since many 
dimers are required to destabilize crystal structures, in 
practice, these mixtures may remain liquid-like, phase 
separate, or form other ordered structures not considered 
here. Along this same line, it is possible that, in specific 
cases, systems may form DQCs or other FK structures 
based on mobile surface particles alone; the entropic ef- 
fect may be stronger for terminal groups that are longer 
or more complex than the one-bead model tested here; 
however, it seems likely that the A15 structure would 
still demonstrate the strongest entropic response due to 
the minimal surface area mechanism [T^HIl]- Since as- 
phericity is common in many micellar systems that also 
have soft coronas, such as the previously discussed TNS 
micelles, it may not be possible to completely separate 
these two aspects. Our results suggest that even moder- 
ate levels of asphericity may enhance the relative stability 
and/or range of stability of DQCs and approximants for 
systems with squishy surface coatings. 

CONCLUSIONS 

Our results demonstrate a two-part, experimentally- 
feasible assembly strategy for forming 3d DQCs and their 
approximants that can potentially be realized for a wide 
variety of systems. We have introduced three new mod- 
els that form DQCs and/or approximants, including a 
simplified model of a spherical micelle and two tethered 
nanoparticle models that resemble micelle-forming sys- 
tems of dendrimers [HIS] and block copolymers [51151 146j. 
Our study lends strong numerical evidence in support of 



the explanation for the stability of the A15 structure in 
systems of dendrimer micelles [12j [13] and its subsequent 
adoption to help explain the formation of the spherical 
dendrimer DQC [51[T31[TB]. Our results imply that shape 
polydispersity, in addition to surface particle mobility, 
is likely to play a role in stabilizing DQCs and approxi- 
mants in micellar systems observed in experiment. In the 
future, our assembly strategy may be employed to facil- 
itate the design of new systems that can form DQCs at 
the nano- and micro-scale, including dendrimers [^[Sllllj. 
surfactants, block copolymers [51121113, and core-satellite 
nanoparticles [lTtii5] . Our results also suggest that mix- 
tures of spheres and dimers |57H59j might, even without 
surface particle mobility, stabilize DQCs or approximants 
under certain conditions, possibly providing a trivial de- 
sign rule for forming these structures. In addition to the 
implications regarding DQC assembly, our results illus- 
trate a powerful design approach for assembling struc- 
tures by controlling particle shape and functionality to 
mimic the key features of pair potentials ^20) . This paves 
the way for future studies based on mapping complex in- 
teraction potentials to packing models, which can poten- 
tially render currently unrealizable systems experimen- 
tally feasible, or expand the breadth of unique structures 
to more general classes of systems. 

MATERIALS 
Simulation of Model Spherical Micelles 

We perform molecular dynamics simulations in the 
canonical ensemble (constant number of particles, vol- 
ume and temperature) of the model spherical micelles 
(MSMs) using the LAMMPS molecular simulation pack- 
age [38] with periodic boundary conditions and reduced 
Lennard- Jones (LJ) units [53]. The Nose-Hoover ther- 
mostat is used with T=1.0 and timestep = 0.005. The 
simulations that we report in this work are performed at a 
nominal volume fraction, (j) « 0.54, computed assuming 
micelles are space-excluding spheres of diameter 5.27tT. 
(j) « 0.54 is chosen because it is comparable to the volume 
fraction of statepoints observed to form spherical micelles 
in previous simulations of tethered nanospheres [22j. A 
cursory sensitivity analysis, performed by starting from 
an assembled sigma phase at fc=5 and either slowing in- 
creasing or slowly decreasing the box size predicts that 
the sigma phase is physically stable within a range of 
nominal volume fractions 0.47 < (f> <0.57, for the spe- 
cific MSM parameters used in this study. In practice, 
the effective micelle volume fraction will be lower than 
the nonimal volume fraction, since the micelles have a 
bumpy, soft corona and will overlap to form dimers. For 
example, if we account for dimers, the volume fraction is 
reduced from 0.54 to 0.46 for a system with /dimor = 0.25, 
assuming dimers are sphero-cylinders of length 6.72 a 
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and diameter 5.27a. In Fig. 5 we provide a schematic 
of the model spherical micelle (MSM), depicting the ap- 
proximate micelle corona and core sizes, as well as the 
dimer to monomer size ratios, for the MSM simulations 
investigated in this paper. 

As discussed in the main text, surface particles inter- 
act via the purely-repulsive Weeks-Chandler-Anderson 
(WCA) potential, meant to capture excluded volume. 
The WCA potential follows the form gO]: 



U 



4e 



(r - a) < Tc 



_ (r — (r — a)^ ^ 

,0, (r-a)>rc 

(1) 

where Us = — e and the interaction cutoff Tc, is set to 
2^/^. a is an adjustable shifting parameter, set to zero 
here. Surface particles are held to their scaffold sites 
with harmonic springs, defined as U{r) ~ kr"^ , following 
the convention in LAMMPS [3H]- The general simula- 
tion procedure is as follows: simulations are initialized 
by creating a random arrangement of MSMs under di- 
lute conditions (0 < 0.1) with fc=10; the system is then 
slowly compressed until the target box size is reached. 
Starting from the target box size, MSMs are simulated 
with k=2 and further allowed to disorder, k is then incre- 
mentally increased by 0.25 until the final value is reached 
(typically, fc = 4 to 5). For each value of fc, the system 
is run for 10 to 50 million timesteps (large systems are 
run for longer than small systems). Simulations are typi- 
cally run for between 50 and 500 million total timesteps. 
This procedure of slowly increasing k mimics the pro- 
cedure used to simulate tethered nanospheres, discussed 
below. Note, our simulations did not span long enough 
time scales to observe tile rearrangements in the ordered 
solid phases. 



Simulation of Tethered Nanospheres 

We perform Brownian dynamics simulations in the 
canonical ensemble of the TNS system with periodic 
boundary conditions and reduced LJ units. The mono- 
TNS building block used consists of a chain of 8 spheri- 
cal beads of diameter a, connected via finitely extensible 
non- linear elastic (FENE) springs pTl 122) . Each chain is 
attached with a FENE spring to a larger "nanosphere" 
of diameter D — 2.5a. The potential energy of the FENE 
spring is given by 



FENE 



(r) 



1 



kRl In 



a) 



Ro 



(2) 



where fc is the spring constant, r is the interparticle sepa- 
ration, Ro is the maximum allowable separation, and a is 
an adjustable shifting parameter to account for excluded 
volume of the nanosphere. For this study, fc=30 and Ro— 




FIG. 5: (a) Histogram of the radial distance between the 
MSM centroid and its mobile surface spheres at fc=5. (b) 
Based on (a), a schematic of the relative size ratio of the 
corona to micelle core in the MSM, where we find the mi- 
celle core diameter is roughly 4cr and the corona extends to a 
diameter of 7cr. (c) The radial distribution function, g(r), cal- 
culated between the centroids of MSMs arranged in a sigma 
structure at fc=5. Note the first peak at ~2.75cr corresponds 
to the center-to-center distance of two overlapping micelles 
(i.e., a dimer). The second peak at ~6(t corresponds to the 
spacing between non-overlapping micelles, (d) Schematic of 
the monomer and dimer size ratios, as derived from g(r) in 
(c). All systems at nominal volume fraction = 0.54. 



1.5, and a — 0.75 for the bond connecting the chain to the 
nanosphere, and zero otherwise. Tethers are treated as 
"solvent-phobic" and thus aggregate at sufficiently low T. 
To model this aggregation, the attractive LJ potential is 
used, give by Eqn. [T] but with Us set to the energy at the 
cutoff and Tc = 2.5. All other interactions are modeled 
by the purely repulsive WCA potential (Eqn. [T] with 
Us — — e and = 2^/^), appropriately radially-shifted 
to account for excluded volume; for tether-nanopshere 
interactions a = 0.75cr and for nanosphere-nanosphere 
interactions a = 1.5cr. Simulations are performed using 
the Brownian dynamics thermostat, where the volume 
fraction of the excluded volume of the individual beads 
is varied between 0.25 < (f> < 0.30, the range where spher- 
ical micelles were predicted in previous work !22]. This 
translates to a nominal micelle volume fraction of ~0.52, 
calculated assuming a characteristic diameter of I2a for 
the spherical micelles (approximated from the radial dis- 
tribution function for micelle centers). 

Di-TNS are modeled in much the same way as mono- 
TNS described above. Chains composed of 4 beads of 
diameter a are connected via FENE springs (Eqn. [2|. 
Two chains are connected to a single nanosphere of di- 
ameter D — 2.0(7, diametrically opposed. This planar 
angle of 180 degrees between the chains is maintained 
by the use of a harmonic spring between the first beads 
of the two polymers, with fc=30 and equilibrium separa- 
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asphericity, as 

FIG. 6: Asphericity histogram of the sigma phase formed 
by the di-tethered nanosphere system. Two representative 
micelles are inset in the plot along with an image of the sigma 
phase formed by the di-tethered nanospheres. 



tion set to Scr. The two polymer chains are chemically 
distinct. One chain is considered to be solvent-phobic 
(i.e., attractive), and thus is modeled by the LJ poten- 
tial. The other chain is considered to be solvent-philic 
(i.e., non-attractive) and modeled by the WCA poten- 
tial. Nanosphere-nanosphere interactions are modeled 
with the LJ potential, appropriately radially-shifted to 
account for excluded volume {a — 1.0). All other in- 
teractions are modeled by the WCA potential, appro- 
priately radially-shifted (a = 0.5 for tether-nanosphere 
interactions). Simulations are performed at = 0.20, 
as calculated from the excluded volume of the individual 
beads. In Fig. 6 we plot the asphericity histogram of the 
micelles formed by the di-tethered nanosphere system. 
We note the asphericity histogram closesly matches the 
result obtained for the mono-TNS system shown in Fig. 
3d of the main text. 

The general simulation procedure used is identical for 
both mono- and di-TNS. Systems start from a disordered 
mixture of TNS, well above the order-disorder tempera- 
ture where little-to-no aggregation occurs (T=2.0). Sys- 
tems are then incrementally cooled to their final tem- 
perature (Twl.O), where systems are run for several mil- 
lion timesteps at each incremental temperature. The po- 
tential energy is monitored to ensure a steady state is 
reached before additional cooling. As the temperature is 
reduced, individual TNS slowly aggregate into micelles; 
micelles form ordered structures at sufficiently low T. 
Simulations are typically run for a total of 40 million 
timesteps. Multiple independent cooling sequences are 
performed to ensure reproducibility of results. Simula- 
tions of mono-TNS are performed for systems of 2500 
building blocks in cubic boxes (22500 total beads) and 
5000 building blocks in boxes with aspect ratio 2:2:1 
(45000 total beads). Simulations of Di-TNS are per- 
formed in cubic boxes for systems of 2000 building blocks. 
See References piHSl] for additional details regarding the 



simulation of TNS. Note, our simulations did not span 
long enough time scales to observe tile rearrangements 
in the ordered solid phases. 

Free Energy of Model Spherical Micelles 

We investigate how surface particle mobility affects the 
stability of MSMs in various crystal structures (Fig. 4a) 
by calculating the change in Helmholtz free energy F 
as a function of the strength of the harmonic springs 
(i.e., k) that anchor particles to the micelle surface, using 
free energy perturbation (FEP) [15]. For each structure, 
the calculation is split into eight independent stages to 
avoid asymmetric bias [50j. consisting of equilibrium sim- 
ulations with spring constants k — 3, 4. ..10. Within each 
stage, we compute the ensemble average: 

F{k,)-Fih) = ~kBT\n(exp (- ^';~/'' ))^ ' (3) 

which gives the free energy difference between stage i and 
j, where j = i ± 1. Each stage consists of a molecular 
dynamics simulation |38j in the canonical ensemble at 
T = 1.0 and micelle nominal volume fraction « 0.54 
to match the simulations. The vertical offset between 
the curves for different structures is computed using a 
FEP variant of the standard Frenkel-Ladd lattice cou- 
pling expansion thermodynamic integration (TI) method 
for molecular systems [52l [53]. We use FEP to adapt 
the method to the complex objects considered here. Al- 
though this method is non-standard, it gives reasonable 
estimates of F that are consistent with our self-assembly 
results. We note this calculation only effects the vertical 
offset of the curves and not how F changes as a function 
of A:. 



Free Energy of Sphere/Dimer Mixtures 

The free energy for mixtures of WCA spheres and 
dimers is evaluated using a three-step scheme based on 
the standard Einstein crystal TI method for spherical 
particles [53H55] . plus an additional alchemical [5^ step 
to convert a given fraction of the spheres into dimers. 
Since systems of WCA spheres do not act like harmonic 
crystals for many of the structures tested, we use the 
Dzugutov (DZ) [T7] system as a convenient reference sys- 
tem that gives harmonic behavior. Computing the work 
required to change the particle interactions from the DZ 
potential to the WCA potential constitutes the third step 
of our scheme. 

In the first step, we compute the free energy dif- 
ference between a non-interacting harmonic (Einstein) 
crystal and a system of spherical particles interacting 
with the DZ potential using the standard Frenkel-Ladd 
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method fS3'-'55'. We denote this free energy difference 

In the second step, we compute the work required to 
change a given fraction, of the spherical particles in 
the system into dimers. We consider a system with the 
energy function: 



t/(A)n = (1 - A)[/p,,e + AC/„ 



(4) 



The free energy required to change the system of DZ 
spheres (pure) to a mixture of DZ spheres and dimers 
(mix) is the integral over the derivative with respect to 
the so-called switching parameter A: 



6X L, 



purc/;^ 



(5) 



The third step is to compute the work required to change 
the potential from DZ to the WCA potential. We con- 
sider a system with the energy function: 



C/(A) 



III 



(1 - A)C/dz + A[/wcA- 



(6) 



The free energy difference for changing the interaction 
potential is given by integrating over the derivative with 
respect to the switching parameter A: 



(7) 



For all calculations, we run 20 independent MC simula- 
tions for different values of A to estimate (5J7(A)/(5A, and 
obtain AF by numerical integration. Simulations are se- 
lectively carried out where |<5^[/(A)/(5A^| is the largest. 
For systems that contain dimers, each simulation begins 
with a compression run with particles constrained to their 
lattice positions allowing rotations and swaps before equi- 
librating at constant density. The total free energy for a 
sphere/dimer mixture is given by: 



F = FEin + AFi + AFii + AF] 



HI- 



(8) 



This formula is used to evaluate F for the WCA system, 
shown in Fig. 4b. 



Isosurface Generation 

Representative configurations from self-assembly simu- 
lations are plotted as time-averaged isosurfaces to coarse- 
grain over thermal fluctuations and produce a clearer pic- 
ture of the structure. To generate the isosurfaces, we 
replace the centroids of the MSMs (or aggregating teth- 
ers in the case of TNS) with a normalized Gaussian of 
width 1.5(7, mapped to a voxel grid composed of cells of 
length Icr, to achieve a degree of spatial coarse-graining. 
We then average the voxel data for ten configurations 
generated within a time window that is much shorter 



than diffusion timescales. This voxel data is then used to 
create isosurfaces within the Visual Molecular Dynamics 
software program [48] , with an isovalue typically ranging 
between 0.1 and 0.2 for MSMs and between 4 and 6 for 
the TNS systems. 
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